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The  Mean  and  Variance  of  Diffuse  Scattered  Power 
as  a  Function  of  Clutter  Resolution  Cell  Size 


1.  INTRODUCTION 


In  general,  analyses  of  the  statistical  properties  of  EM  signals  scattered  from  random 
rough  surfaces  have  assumed  that  the  surface  contributing  to  the  scattering  Is  infinite 
compared  to  the  characteristic  dimension  of  the  roughness  (In  particular,  the  correlation 
loigth  T).  Further,  the  surface  Is  assumed  to  be  homogeneous  and  Isotropic.  The  result  of 
these  assumptions  is  that  the  scattered  signal  tends  to  have  a  Rayleigh  distribution.  The  mean 
scattered  power  Is  the  product  of  the  mean  power  per  unit  area  and  the  area  of  the  cell.  If  the 
area  Is  reduced  by  shrinking  the  range  extent  of  the  ceO  (corresponding  to  shorter  signal 
pulses)  then  the  mean  clutter  power  scattered  from  the  surfiice  Is  reduced  so  long  as  the 
statistical  distribution  remains  unchanged.  This  report  addresses  this  phenomenon  by 
examining  the  effect  of  finite  cell  size  on  the  calcvilated  statistical  properties.  The  approach  Is 
to  analyze  the  behavior  of  the  nomudlzed  mean  cross  section  o**  and  normalized  variance  o*^ 
of  the  power  as  a  function  of  cell  size.  In  this  Initial  stu^  the  surfoce  roughness  Is  assumed  to 
be  one-dimensional. 

First,  the  theoretical  model  will  be  '^escribed,  followed  by  an  outline  of  the  numerical 
techniques  used  to  derive  the  results.  Next,  the  main  part  of  the  report  presents  the  results  and 
analyzes  them.  The  behavior  of  o**  and  o*”  Is  shown  as  a  function  of  scattering  angle  0,.  Plots 
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of  0“  and  a"*  vs.  0^  are  studied  for  selected  values  of  ai^e  of  Incidence  0j ,  clutter  cell  size  L, 
rms  surface  slope  V2  o/T  and  EM  wavelength  (carrier  wavelength).  The  two  cases  where  the 
surface  properties  are  either  a  perfect  conductor  or  a  los^r  dielectric  are  Investigated.  The 
plots  of  a°°  vs.  0^  compare  the  actual  value  of  o“*  to  the  value  of  0“°  assuming  the  scattered 

power  is  a  Rayleigh  distribution.  In  the  final  section,  the  results  are  summarized  and 
conclusions  ofiered. 


2.  THEORETICAL  FORBOULATION 


The  analysis  begins  by  assuming  that  the  frequency  of  the  EM  wave  is  characterized  by  the 
carrier  frequency  co.  and  that  the  corresponding  wavelength  X  ■  3  K£  is  much  smaller  than  the 

CD 

rough  surface  correlation  length  T,  that  is.  T  »  X.  Then,  the  wave  scattering  process  may  be 
described  by  physical  optics  (PO).  Both  the  correlation  length  and  the  surface  heights  are 
assumed  to  be  described  by  Gaussian  functions.  The  rough  surface  heights  vary  in  one 
dimension  only,  4  *  z  (x).  Then,  following  the  development  given  in  a  previous  report*,  the 
scattering  takes  place  In  the  plane  of  variation  (^^  »  Of*)  and  there  is  no  shadowing. 

The  normalized  cross  section  of  a  rough  surface  o®  Is  proportional  to 

<,”  .  <  (E./E,)  .  (£./£,)•>  -  <(£,/£,).'  II, 


where  <  •  >  denotes  an  ensemble  average  over  the  variables  ^j,  Pj,  Jij-  Here,  E^  is  the 
scattered  electric  field,  Ej  is  the  incident  electric  field,  wid  ^  are  the  random  surface 
heights  at  two  points  Xj  and  Xj  and  pj  and  random  surfoce  slopes  at  the  two  points. 

The  most  general  expression  for  o®  for  a  one-dimensional  rough  surface  involves  a  six-fold 
integral  over  the  variables  Xj.  x^,  p^.and  p^.  By  using  the  fact  that  the  surface  heights 

are  to  be  regarded  as  a  stationary  random  process,  so  that  the  correlation  function  of  the 
surface  heights,  c  (t).  Is  a  function  only  of  the  separation  between  points  to  Xj  -  Xj  ,  the  six¬ 
fold  integral  for  0°  may  be  reduced  to  a  five-fold  integral.  Then,  using  the  fact  that  the  general 
expression  for  o®  is  a  function  only  of  the  hel^t  differences  §  »  ^  .  the  five-fold  integral 

can  be  reduced  to  a  four-fold  Integral.  Making  the  additional  assumption  that  the  rms  slope 
V2  o/T  <  1  (where  o  ■  rms  surface  height)  permits  the  surface  heights  to  be  decorrelated  from 
the  surface  slopes.  Then,  the  integrations  over  the  surface  slopes  pj  and  Pj,  may  be  readily 

performed,  so  that  the  expression  for  o®  for  a  finite  clutter  cell  reduces  to  a  double  Integration. 
Under  these  same  assumptions,  and  using  the  expression  for  the  scattered  electric  field  based 
upon  physical  optics  (PO).  an  expression  for  the  normalized  variance  of  the  power  o®®  may  be 

*  Papa.  R.J.,  and  Woodworth.  M.  (1988)  The  Numerical  EXxiluafion  of  a  Physics  Optics 
Normalized  Cross  Section  for  a  Rough  Surface,  RADC-TR-87-280,  ADA198919. 
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derived.  For  a  clutter  cell  of  finite  size,  the  normalized  variance  a°°  may  be  expressed  as  a 
four-fold  Integral. 

Since  the  surface  Is  assumed  to  be  rough  In  one  dimension  only,  the  scattered  field  Is  a 
cylindrical  wave 


E^  =  V2/  R')  exp[lkj«  R' 


(2) 


where  =  scattered  wave  number 

R'  =  R„-r 

Rjj  =  distance  from  origin  to  field  point 
r  =  distance  fi-om  origin  to  source  on  the  surface. 

The  scattered  field  is  calculated  evaluating  the  Helmholtz  Integral  (Beckmann  and 
Splzzlchlno^): 


where 

V  =  V2/ (nk^  •  R')  exp[lk^«R']  (4) 


—  »  normal  derivative  on  the  surface  S.  (Here,  we  Integrate  over  only  one  dimension,  x). 


^  Beckmann.  P..  and  Splzzlchlno,  A.  (1987)  The  Scattering  of  Electromagnetic  Waves 
from  Itough  Surfaces.  Artech  House,  Inc..  Dedham,  MA. 
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It  is  assumed,  under  the  PO  approximation,  that  the  totld  field  on  S  is  given  by 


E  =  (1  +  R)  E, 


(5) 


where 


R  =  Fresnel  plane  wave  reflection  coefficient. 


By  combining  Eqs.  (2),  (3).  (4).  and  (5)  and  using  the  assumptions  given  in  the  previous 
discussions,  an  expression  for  the  scattered  wave  may  be  obtained  (see  Beckmann  and 
Splzzlchlno^): 


where 


V  »  k  [(sin  9|  -  sin  9^)  x 


(cos  9j  +  cos  9^)  z] 


=  V  X  +  V  z 
X  z 


and  k  =  2  k/\. 


Here. 


F  » [2R  (l  +  cos  (9j  +  9J)  /  (cos  9^  +  cos  9j] 


(7) 


and  we  have  assumed  a  horizontally  polarized  incident  wave,  and  thus  we  also  specify  the 
appropriate  form  for  R. 
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The  expression  for  the  normalized  cross  section  a°  of  a  one-dlmenslonal  rough  surface  Is 
given  by; 


.  (2  n  R,/L)  I  <  iE,/EjS  -  <  |E./E,!  >’] 


(8) 


where  <  •  >  denotes  an  ensemble  average  over  the  surface  heights.  One  may  then  derive  the 
following  expression  for  o°: 


O 

a  = 


2nLX) 


ru/2 

|.IV2 

dx, 

-L/2 

f-U2 

I  cos  v^  t 


where  x  =  Xj  -  x^ 


(9) 


Xi 


characteristic  function  for 


exp 

characteristic  function  for 


bivariate  height  distribution 
univariate  height  distribution 


Rayleigh  parameter  squared. 


The  normalized  variance  In  clutter  power  Is  defined  In  terms  of  the  fourth  moment  of 
the  scattered  field: 


OO 


a 


|e./e/  >“). 


(10) 


First,  consider  the  term 


•U2  -L/2 

<!EyEj^>oc  dXj  dx2e‘’^«<*»“*2> 

J-L/2  J-L/2 


(11) 
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d*ae*PK(*i-*a)]X2- 


Next,  consider  the  square  of  the  second  moment  of  the  scattered  wave: 


<|E,/Ef  ><|E,/EJ^  > 


(12) 


.L/2  aU2  0U2  ^L/2 

“  dxj  dxj  dx3  <1x4  •  «P[iV,((Xj-Xj)  +  (X3-xJ)] 

J-L/2  J-U2  J-l/2  J-l/2 


The  ensemble  average  of  the  fourth  moment  of  the  scattered  dectrlc  field  is  given 


fVa 

pl/2 

pU2 

/>«. 

dXj 

1  dx. 

dXg 

% 

f-l/2  « 

f-U2  J 

(-L/2  « 

'-L/2 

l-Xg-X^) 


(13) 


Define  the  fourth  order  characteristic  function  X4  *1* 


X  ^  -  <  e‘^* > 


(14) 
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=  j  '*«2|  <*«3j 


where  P  ^4)  fourth  order  height  distribution  function.  For  a  Gaussian 

distribution. 


P  (4,.  ^2.  y  =  ((2n)*[DET  C]‘^)  *  exp  [-(1/2)  ^  C"^  ^ 


(15) 


where 


^3 


L  -^4 


and  the  covariance  matrix  C  Is  given  by 


1 

C12 

C.3 

Cu 

II 

C21 

1 

C24 

C3, 

^32 

1 

C34 

1^41 

C42 

C43 

1 

It  Is  known  that  the  Fourier  transform  of  a  Gaussian  Is  also  a  Gaussian  (Stark  and 
Woods^): 


I  exp 


-(l/2)a)^C  (0 


(16) 


^  Stark.  H.,  and  Woods,  J.  (1986)  ProbcdfUtty,  Random  Processes  and  Estimation  Theory 
for  Engineers,  Prentice-Hall.  New  York. 
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where 


After  some  algebraic  manipulations,  the  expression  for  X4  may  be  derived: 

*4  *  exp  [-^  (2  -  C12  +  Cj3  -  -  C33  +  -  C^)].  (17) 

Ihen,  combining  Eqs.  (lOl  through  (16),  one  may  obtain  the  following  eiqjresslon  for  the 
variance  In  power 

2  eL/2  fU2  Ahn  aU2 

o~.(f^/[2xXL])  dx,  dXj  dXg  dX4«)S[v^(X,-X3  +  X3-X^)]*[X4-X,2X34].(18) 
J-L/2  J-IJ2  J-U2  J-l/2 

Here.  X12  =  op  (1  -  CjjM 
and  X34  *  exp  I- 1*  (1  -  Cg^)). 

If  the  statistical  distribution  of  the  power  scattered  imm  a  clutter  cell  is  Haylel^.  then 

o®«P,  (19) 


vdiere  Is  the  mean  of  the  scattered  wave.  Eq.  (19)  Is  true  for  any  statistical  distribution.  But 
the  variance  In  the  power,  for  a  Rayleigh  distribution.  Is  given  \jy\ 


(20) 


Eq.  (20)  Is  true  otify  for  a  Raleigh  distribution.  The  deftnltlons  of  0°  and  o°°  given  In  this 
report  are  such  that,  for  a  Raykfgh  distribution: 
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3.  NUMERICAL  METHODS 


The  four-fold  Integral  in  Eq.  (18)  and  the  two-fold  Integral  In  Eq.  (9)  must  be  solved 
numerically.  Any  quadruple  Integral  will  require  many  hmctlon  evaluations.  The  key  Is  to 
use  an  efficient  numerical  Integration  algorithm  that  minimizes  the  required  number  of 
function  evaluations.  We  used  three  different  methods:  two  Monte  Carlo  routines  and  an 
adaptive  Gauss  quadrature  routine.  These  three  methods  have  different  strengths  and 
weaknesses,  thus  none  Is  superior  to  the  others  In  all  situations. 

The  first  Monte  Carlo  routine  we  used  Is  a  modified  version  of  a  routine  obtained  from  the 
Press  and  Farrar^  paper,  titled  "Recursive  Startlfied  Sampling  for  Multidimensional  Monte 
Carlo  Integration".  The  four-dlmenslonal  region  to  be  Integrated  Is  first  sampled  to  determine 
the  distribution  of  the  function.  The  region  Is  then  divided  up  Into  subregions  each  assigned  a 
certain  percentage  of  the  points  where  the  function  is  evaluated.  In  this  way.  the  laigest 
number  of  function  evaluations  are  concentrated  where  the  function  is  most  significant.  The 
points  In  each  four-dlmenslonal  subregion  are  determined  fay  a  uniform  random  number 
generator.  The  function  Is  evaluated  at  these  points  and  the  Integrand  value  for  each 
subregion  is  determined.  The  overall  Integral  Is  found  combining  all  subregions.  The 
standard  deviation  of  the  answer  is  also  calculated.  If  the  standard  deviation  is  too  large 
(>  1%),  then  the  number  of  points  is  increased  up  to  a  maximum  of  about  7  x  10^.  B^nd  this 
number  the  round-off  error  becomes  the  main  factor  limiting  the  accuracy. 

The  second  Monte  Carlo  routine  Is  based  on  the  Press  and  Teukolsl^  paper  titled  "Quasl- 
(that  Is.  Sub-)  Random  Numbers".  Simple  Monte  Carlo  Integration  relies  on  points  distributed 
at  random  over  the  region  to  be  integrated.  The  fitictlonal  error  of  the  estimated  integral 
decreases  1/Vn.  where  N  is  the  number  of  sample  points.  Press  and  Teukolsl^  show  that 
when  the  multidimensional  points  are  determined  a  subrandom  Sobol  sequence,  the 
fractional  error  decreases  at  a  rate  close  to  1/N.  (The  points  generated  the  Sobol  sequence 
appear  random  when  plotted,  but  actually  the  Sobol  sequence  deterministically  chooses  points 
In  space  which  maximally  avoid  each  other  and  fill  the  space  more  uniformly  than 
uncorrelated  random  points.)  The  Integral  can  be  accurate^  determined  with  mai^  fewer 


*  Press.  W.H..  and  Farrar.  G.R.  (1990)  Recursive  stratified  sampling  for  multldlmenslcmal 
Monte  Carlo  Integration.  Computers  (n  Physics,  4:190-195,  March/Aprll. 

^  Press,  W.H.,  and  Teukolshy,  S.A.  (1989)  Quasi-  (that  Is,  sub-)  random  numbers.  Computers 
In  Physics.  3:76-79,  November/December. 


function  evaluations  than  the  standard  Monte  Carlo  method.  The  function  is  evaluated  at 
each  point  and  the  values  are  summed  to  give  the  value  of  the  Integral.  The  standard  deviation 
Is  also  detennlned.  If  the  result  is  not  accurate  enough,  the  Sobol  sequence  Is  simply  extended 
until  the  desired  accuracy  is  reached,  up  to  a  limit  of  2^*  points. 

The  adaptive  Gauss  quadrature  routine  is  based  on  an  Int^atlon  routine  presented  by 
Shamplne  and  Allen^  We  converted  their  routine  from  one  based  on  Simpson's  rule  to  one 
based  on  four-point  Gauss  quadrahme.  We  further  modified  the  routine  to  solve  the  multi¬ 
dimensional  integration,  and  to  Improve  its  efilclency.  The  basic  technique  is  simple.  First 
an  integral  is  evaluated  with  four-point  Gauss  quadrature,  then  the  integral  is  bisected  and 
each  half  Is  evaluated  separately  with  four-point  Gauss  quadrature  and  the  halves  summed.  If 
the  two  results  are  different  by  more  than  the  desired  accuracy,  then  the  right  half  is  stored 
and  the  left  half  Is  again  bisected  and  each  new  section  Is  evaluated.  Again  the  accuracy  Is 
checked  by  comparing  the  sum  of  the  evaluations  over  the  two  quarters  with  that  of  the  half 
Interval.  Bisection  continues  until  either  the  desired  accuracy  Is  reached  or  the  subsection 
becomes  too  small.  In  either  case,  the  values  of  the  subsections  are  stored  and  the  routine 
moves  to  the  next  subsection  to  the  right,  adiich  Is  also  bisected  until  the  desired  accuracy  Is 
reached.  In  order  to  maintain  the  given  overall  accuracy,  the  accuracy  Is  divided  by  the  square 
root  of  2  with  each  finer  level  of  bisection.  The  evaluations  over  the  set  of  subsections  are 
summed  to  give  the  total  integrand.  This  adaptive  method  concentrates  the  computational 
effort  vdiere  It  is  most  needed  and  avoids  wasting  effort  where  It  is  not  required.  This  Is  done 
for  all  four  levels  of  Integration.  In  the  outer  Integration,  an  error  threshold  Is  detennlned 
which  Is  passed  to  the  inner  integrals.  If  the  value  of  a  subsection  falls  below  the  threshold,  it 
Is  marked  "done"  even  if  It  has  not  converged,  because  the  value  Is  too  small  to  significantly 
contribute  to  the  final  result. 

The  three  methods  are  driven  by  different  mechanisms.  The  number  of  function 
evaluations  is  a  user  input  for  our  first  Monte  Cario  method  and  the  enor  bound  of  the  result 
Is  determined  from  the  data.  Since  the  conq>uter  CPU  time  is  proportional  to  the  number  of 
evaluations,  we  know  ahead  of  time  how  long  the  first  Monte  Carlo  method  will  take. 
Conversely,  the  error  bounds  are  a  user  Input  for  the  second  Monte  Carlo  method  and  the 
adaptive  Gauss  quadrature,  and  the  function  is  evaluated  as  many  times  as  necessary  to 
achieve  that  accuracy;  thus  the  amount  of  time  can  vary  widely  depending  on  the  exact 
structure  of  the  function.  The  first  Monte  Carlo  method  Is  also  susceptible  to  accuracy 
problems  If  the  function  oscillates.  The  adaptive  methods  are  better  able  to  compensate  for 
any  oscillations.  The  oscillations  come  primanly  from  the  term  cos  (v^  (Xj  -  X2  X3  -  x^))  in 

Eq.  (18).  Doubling  the  length  of  a  cell  (L)  will  double  the  range  of  each  dimemslon  of  the 
integral  and  thus  double  the  number  of  oscillations  in  each  dimension.  Another  factor  is 
Vjj  =  k  (sin  0,  -  sin  which  contributes  the  effect  of  wavelength  and  the  angle  terms.  The 
further  apart  the  Incident  and  scatter  angles,  the  more  oscillations.  Generally,  when  v^  is 
small,  any  method  Is  accurate,  but  the  Monte  Carlo  methods  are  much  more  efficient,  up  to 
orders  of  magnitude  faster  than  adaptive  Gauss  quadrature.  As  the  number  of  oscillations 


®  Shamplne,  L.F..  and  Allen.  R.C.  (1973)  Numerical  Confuting:  An  Introduction,  W.B. 
Saunders  Company. 
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Increase,  the  Hrst  Monte  Carlo  method  becomes  Inaccurate  and  the  second  Monte  Carlo 
method  and  Gauss  quadrature  become  the  more  accuracte  methods  of  integration.  However, 
the  run  time  of  the  Gauss  quadrature  routine  also  Increases  dramatical^  as  the  oscillations 
Increase.  These  two  factors  put  a  practical  limit  on  the  range  of  scatter  angles  that  can  be  used 
with  a  given  data  set.  The  second  Monte  Carlo  method  Is  generally  the  fiastest  and  most 
accurate  method  for  the  widest  range  of  angles. 

The  two-fold  Integral  Cq.  (9)  Is  evaluated  In  a  similar  manner  as  the  four-fold  Integral  (Eq. 
(18)].  As  expected,  the  two-fold  Integral  can  be  evaluated  much  faster  than  the  four-fold 
Integral.  Once  the  Integral  Is  found  It  Is  a  simple  matter  to  convert  from  o°  and  o*^  for  a 
perfect  conductor  to  o°  and  a°°  for  a  dielectric  surface.  Slmpty  use  the  appropriate  value  of  F 
in  Eqs.  (9)  and  (18). 


4.  ANALYSIS  OF  NUMERICAL  RESULTS 


In  this  section,  a  number  of  graphs  will  be  presented  of  o'*  vs  6^  and  0°°  vs  6^,  with 
parameters  that  Include  the  angle  of  Incidence  6j,  wavelength  X,  standard  deviation  In  surface 
height  o.  surface  correlation  length  T,  length  of  rough  surface  L,  and  complex  dielectric 
constant  e.  The  curve  Indicated  with  long  dashed  lines  represents  o**  as  calculated  numerically 
according  to  Eq.  (9).  The  curve  Indicated  by  crosses  represents  o”®  as  calculated  numerically 
according  to  Eq.  (18).  The  curve  Indicated  ty  short  dashes  represents  o“®  for  a  Rayleigh 
distribution,  as  given  by  Eq.  (21)  with  the  value  of  o®  given  Eq.  (19). 

In  Figures  1  through  24.  we  will  study  the  behavior  of  o®  and  o®®  vs  scattering  angle  6,  for 
different  parameter  regimes.  First,  the  angle  of  incidence  6,  Is  chosen  to  be  relatively  small 
(6,  =  30®)  and  then,  later,  is  taken  to  be  large  (6,  >  75®).  The  wavelength  Is  first  taken  at  L-band 
(X  s  0.25m)  and  then  decreased  to  S-band  (X  >  0.1m).  The  cell  size  Is  first  taken  to  be  small 
(L  3  3T),  then  Increased  to  L  »  €T,  and  eventually  to  L  s  12T.  Also,  the  surface  roughness  starts 
out  at  a  relatively  small  value  (o/T  s  0.2)  and  then  Is  increased  to  o/T  *  0.5. 

In  Figure  1 .  the  values  of  the  parameters  are  9, «  30®.  X  *  0.25m,  o  >  0.2,  T  «  Im,  L  =  3m,  and 
e  Is  a  perfect  conductor.  For  a  large  clutter  cell  (L  >  lOT).  the  scattered  power  Is  Ra^el^ 
distributed  (Beckmann  and  Splzzlchlno^).  The  comparison  of  the  exact  o®®  with  the  Rayleigh 
value  of  o®®  shows  that  for  a  short  cell  size  (L  »  3T),  the  clutter  power  Is  Rayleigh  distributed  In 
an  angular  region  about  the  specular  direction  (9,  ■  9|).  where  the  angular  region  extends  fiom 
about  9,  B  50®  to  0®.  As  9^  moves  away  fix>m  the  angular  region,  the  clutter  power  becomes 

more  non-Rayleigh,  and  the  departure  from  the  Rayleigh  distribution  Is  more  dramatic  In  the 
backscatter  region  9^  <  0®  than  In  the  forward  scatter  region  9,  >  50® . 

In  Figure  2.  the  parameters  are  the  same  as  In  Figure  1.  except  that  now  e  is  a  complex 
dielectric  (30  -i-  J2).  representing  moist  loam.  Although  the  results  will  not  be  shown  here,  a 
whole  set  of  complex  dielectric  constants  were  Investigated,  where  the  parameters  were  the 
same  as  In  Figure  1  except  that  e  ■  (2  +  J4),  (2  +  JO).  (4  +  j2).  (4  +  JO),  (80  +  j9).  and  (80  +  JO).  The 
results  may  be  summarized  by  stating  that  the  values  of  o®  and  o®®  are  less  for  a  complex 
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dielectric  surface  compared  to  a  perfect  conductor.  Indicating  a  general  trend.  Also,  It  may  be 
stated  that,  for  a  complex  dielectric  surface.  Increasing  either  the  real  or  imaginaiy  part  of  e 
causes  both  o°  and  to  Increase.  In  Flgvire  3.  the  parameters  are  the  same  as  In  Figure  1 , 
except  that  the  cell  size  Is  larger,  with  L  s  6T.  Comparing  Figure  3  with  Figure  1 ,  It  may  be 
noted  that  the  larger  cell  size  results  In  a*"  being  closer  to  the  Rayleigh  distribution.  In  Figure 
4,  the  parameters  are  the  same  as  In  Figure  2.  except  that  L  a  6T.  Again,  the  larger  cell  size 
results  In  being  closer  to  the  Rayleigh  value,  for  the  surface  with  a  complex  dielectric 
constant. 

In  Figures  5  and  6.  the  parametoa  are  the  same  as  In  Figures  1  and  2.  except  that  o  ■  0.5m, 
so  that  o/T  «  0.5  (a  rougher  surface).  Comparing  Figures  5  and  6  with  Figures  1  and  2.  one  may 
note  that  o”  is  more  constant  in  value  (flatter)  as  6^  varies,  indicating  that  the  router  surface 

scatters  more  Isotropically,  both  for  a  perfect  conductor  and  a  los^  dielectric.  Also,  o**^  Is  less 
close  to  the  Rayleigh  value  of  (o**)^  as  the  surface  gets  rougher,  both  for  a  perfect  conductor  and 
a  lossy  dielectric.  In  Figures  7  and  6,  the  parameters  are  the  same  as  In  Figures  5  and  6,  except 
that  the  cell  size  is  larger,  L  b  6T.  The  Increase  in  cdl  size  causes  the  true  variance  to 
decrease,  becoming  closer  to  the  variance  for  a  Rayleigh  distribution. 

In  Figures  9  and  10.  the  an^e  of  incidence  has  been  dian^d  so  that  6,  ■  75**.  with  aU  the 

other  parameters  equal  to  those  in  Figures  1  and  2.  Comparing  Figures  9  and  10  with  Figures 
1  and  2.  one  may  note  that  Increasing  the  angle  of  Incidence  results  In  smaller  values  of  0“  and 
for  both  a  perfect  conductor  and  a  lossy  dielectric.  Again,  It  may  be  noted  that  the 
variance  Is  close  to  a  Rayleigh  distribution  in  an  angular  region  about  the  specular 
direction  6^  ■  6,.  and  deviates  from  It  as  0^  moves  away  from  this  angular  region.  As  in  Figures 

1  and  2.  the  departure  from  the  Rayleigh  distribution  Is  more  dramatic  for  the  backscatter 
region,  6^  <  0°  than  for  the  forward  scatter  region. 

In  Figures  1 1  and  12,  the  parameters  are  the  same  as  In  Figures  9  and  10,  except  the  cell 
size  Is  larger,  L  «  6T.  The  Increase  In  cell  size  results  In  the  variance  o*”  values  becoming 
closer  to  the  Rayleigh  distribution  values,  for  both  a  perfect  conductor  and  a  lossy  dielectric. 

In  Figures  13  and  14,  the  parameters  are  the  same  as  in  Figures  1 1  and  12,  except  that  the  cell 
size  has  been  Increased  further,  to  L  >  12T.  This  further  Increase  results  In  the  ty°°  values 
becoming  almost  equal  to  the  Rayleigh  distribution  values. 

In  Figures  15  and  16,  the  parameters  are  the  same  as  In  Figures  9  and  10,  except  that  the 
standard  deviation  in  surface  height  has  been  increased,  o  >  0.5m,  so  that  the  surface  Is 
rougher  with  o/T  >  0.5.  The  Increase  In  svirface  roughness  causes  the  values  of  o”  and  o*^  to 
increase.  Also,  the  vaues  of  0°**  deviate  more  pronouncedly  from  the  Rayleigh  distribution 
values.  In  Figures  17  and  18,  the  parameters  are  the  same  as  In  Figures  15  and  16.  except  that 
the  cell  size  has  been  increased  to  L  s  6T.  As  shown  previously,  the  increase  In  cdl  size  causes 
the  o°"  values  to  be  closer  to  those  of  a  Rayleigh  distribution. 

In  Figures  19  through  24.  the  EM  wavelength  has  been  dianged  to  X  «  0.1m,  so  that  S*band 
microwave  radiation  is  being  considered,  instead  of  L-band  (X  ■  0.25m).  The  parameters  in 
Figures  19  and  20  are  the  same  as  In  Figures  1  and  2.  except  for  the  EM  wavdoigth. 

Comparing  Figures  19  and  20  with  Figures  1  and  2,  one  may  note  that  the  decrease  In  X  results 
In  a  slight  reduction  in  o®  for  0^  <  0®.  Also,  the  o®®  values  have  Increased,  with  a  greater 

deviation  from  the  Rayleigh  distribution  values. 
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In  Figures  21  and  22.  the  parameters  are  the  same  as  in  Figures  9  and  10,  except  that  the 
EM  wavelength  has  been  reduced  to  X.  »  0.1m.  Comparing  Figures  21  and  22  to  Flgiarcs^9  and 
10.  it  may  be  observed  that  the  decrease  in  X  causes  the  a®  values  to  decrease  for  0,  <  0®.  Also. 

the  o®®  values  have  Increased,  exhibiting  a  greater  deviaUon  from  the  Rayleigh  distribution 
values. 

In  Figures  23  and  24,  the  parameters  are  the  same  as  in  Figures  21  and  22,  except  that  the 
ceU  Size  has  been  increased  to  L  =  6T.  As  at  L-band.  the  Increase  in  ceU  size  causes  the  o®® 
values  to  be  closer  to  the  Rayleigh  distribution  values. 
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SUMMARY  AND  CONCLUSIONS 


Reviewing  Figures  1  through  24,  some  general  statements  about  the  behavior  of  and  o'* 
may  be  made.  First,  for  our  definitions,  &***  *  (o'*)^  for  a  Raleigh  distribution.  The 
normalized  variance  is  always  less  than  ttie  mean  cross  section  o".  The  normalized 
variance  0°'*  and  the  mean  cross  section  o'*  both  decrease  as  the  an^e  of  incidence  6j  Increases, 

with  all  other  parameters  fixed.  The  normalized  variance  o**'*  becomes  dose  to  Raylel^ 
scatter  as  the  cell  size  L  Increases.  In  the  specular  direction.  &***  is  very  dose  to  the  Rayleigh 
value  of  (o")^.  The  mean  normalized  cross  section  decreases  as  the  cell  size  Increases  only 
for  the  backscatter  region  (6,  <  O'*).  Both  o'**  and  o*  are  smaOer  In  magnitude  for  a  los^ 

dldectrlc  than  for  a  perfect  conductor. 

There  are  several  factors  controlling  the  distribution  of  the  scattered  power:  (1)  The 
dlfiuse  power  scattered  In  the  region  inside  of  an  angular  band  surrounding  the  specular 

direction  tends  to  be  Ri^ldgh  for  ceU  sizes  L  >  3T.  (2)  When  the  cell  size  L  <  6T.  the 

distribution  of  scattered  power  Is  non-Ri^lel^  outside  of  an  angular  region  surrounding  the 
specular  direction.  6,  +  SO'*  <6,  <  6j  -  30*.  (3)  For  a  fixed  cdl  size  L,  the  distribution  becomes 
less  Rayleigh  as  the  scattering  angle  0,  moved  away  fit>m  the  angular  band  surrounding  the 

specular  direction,  and  the  departure  from  the  Rayleigh  distribution  is  more  pronounced  In  the 
backscatter  direction  0,  <  0*.  compared  to  the  forward  scatter  direction. 

These  trends  are  true  for  both  a  perfect  conductor  and  a  los^  dielectric.  Also,  these  trends 
are  true  at  both  L-band  and  S-band  fipequendes.  Increasing  the  frequency  from  L>band  to 
S-band  causes  the  distribution  of  the  scattered  power  to  become  more  non-Rayleigh- 
Finally,  It  may  be  noted  that  Increasing  the  surface  roughness  tends  to  make  the 
distribution  of  the  scattered  power  more  non-Raylelgh.  This  Increase  In  surface  roughness 
may  be  due  to  either  an  increase  in  the  standard  deviation  of  surface  height  o  or  a  decrease  in 
the  wavelength  Increasing  the  surface  roughness  tends  to  make  the  dlfiuse  scattered  power 
more  spatially  Isotropic,  in  general.  Also,  the  angular  width  of  the  r^lon  in  0^  about  the 
specular  direction  where  o*  has  significant  magnitude  decreases  as  the  angle  of  incidence  0, 
Increases. 

The  distribution  In  scattered  power  tends  to  be  more  Raleigh  in  and  near  the  specular 
direction,  because  there  are  more  facets  which  can  contribute  to  scattering  in  this  direction. 
There  is  a  Gaussian  distribution  In  the  slopes  of  the  fricets,  and  the  maximum  of  this 
distribution  of  slopes  lies  precisely  In  the  specular  direction.  There  are  relatively  few  fricets 
with  ai^MToprlate  slopes  In  the  backscatter  direction  <  0*).  so  that  the  departure  frtnn  the 

Rayleigh  distribution  Is  large  for  the  scattered  power. 
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OF 

ROMS  LABORATORY 

Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re¬ 
search,  development,  test,  and  technology  transition  in  support  of  Air 
Force  Command,  Control,  Communications  and  Intelligence  (C^I)  activities 
far  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  sui^ort  within 
areas  of  competence  is  provided  to  BSD  Program  Offices  (POs)  and  other 
BSD  elements  to  perform  effective  acquisition  of  C^I  systems,  bi  addition, 
Rome  Laboratorys  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  oth&r  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliabUity/maintainability  and  testability. 


